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ABSTRACT 

We present out results on the cross-correlation of Lyman break galaxies (LBGs) around 
three damped Lya absorbers (DLAs) at Zabs — 3 from deep {fJ-i,AB{sky) ~ 27.6mag arcsec"^) 
UBVI KPNO 4m/M0SAIC images. The large area of the MOSAIC images, 0.31 deg^ or 
~ 65 X 65hjl Mpc co-moving at redshift z = 3, allows us to probe the clustering of LBGs on 
scales up to 20 Mpc co-moving. Our survey covers a total of 1 deg^ and contains ~3,000 LBGs 
with photometric redshifts between 2.8 and 3.5. Using the redshift likelihood distributions with 
mj as a prior, we selected LBGs within a redshift slice of width = 0.15 (corresponding to a^, 
the uncertainty in photometric redshifts) centered on the redshift of the absorbers. Within that 
redshift slice, we find that the DLA-LBG cross-correlation Wdg is Wdg = (1.62ibl.32) x Wgg, where 
Wgg is the LBG auto-correlation. This corresponds to a correlation length of To = 5 it A.5h~^ 
(co-moving) (or ro = 7 ± 6.8/if/ Mpc). The cross-correlation is most significant on scales 
5 — 10 Mpc. Through Monte Carlo simulations, we find that Wdg is significantly greater than 
zero at the > 95% level. In three other redshift slices that do not contain a DLA, we do not 



find any evidence of clustering. A larger sample will enable us to discriminate between Wdg/ 
< 1 or Wdg/wgg > 1, i.e to test whether DLA halos are more or less massive than LBG halos. 

Subject headings: cosmology: observations — galaxies: evolution — galaxies: high-redshift — 
quasars: absorption lines — quasars: individual (APM08279-I-5255, 
PC1233-F4752, J0124-h0044) 
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1. INTRODUCTION 

QSO absorption lines, including damped Ly-a ab- 
sorbers (DLAs), and Lyman break galaxies (LBGs) 
are currently our two major sources of information on 
high redshift galaxies. After more than two decades 
of study, the exact nature and detailed characteristics 
of damped absorbers remain unexplained. Here, we 
seek to constrain the properties of DLA halos using 
LBGs (|Steidel et al.l . Il999i ) as tracers of large scale 
structure. 

DLAs contain the largest rese rvoir of neutral hvdro- 
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anzetta et al.l.ll991 
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gen (Hi) at high redsh ifts (e 
Lanzetta, Wolfe. & Tumshe 

(200 li) . They contain more neutral Hi than all the 
absorbers in the Ly-alpha forest combined. Mor- 



ever, the amount of Hi in DLAs at high redshifts 
corres ponds to the amoun t of Hi in stars today at 
z = : lElhson et al.l (l200ll ) find 0///(z = 3) = 10 
while iBell et al.1 (l2003l ) rneasure n^z = 0) = 10" 
(bo th numbers a re fo r h = 0.65). These facts 
led IWolfe et al.l (Il986t ) to put forward the 'disk 
hypothesis', namely that DLAs are large thick 
gaseous disk galaxies. Despite the numerous ob- 
servations directed at DLAs in the past decade 
' e.g. imaging studies such as iMoller fc Warren . 



199,1 : iLg wentha l et al.l . 1199,4 ISteidel et all 
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ibo Moller. WarrenL Il999l: iKulkarni et al. . 
2000"; 'Pettini et al.', '2000'; 'Rao Rr. Turnsheki l200r : 
Bouche et al.. 2001: Meille r et al., 2002), this hypoth- 
esis has been debated and the role of DLAs in galaxy 
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formation is still not unders tood (see the more ex- 
haustive summary of |Pettinil .|20Q3). 
Hints to the natu r e of DLAs are given by numerical 

(iKatz et alJ. Il996l: Haehnelt. Steinmetz. &: Ranch L 

2^0; Oar dner et a1.l .l2nni: N^agamine. Springe l. fc Hernau 
2nn3bl'l and semi-analyt i cal (iKauffniann I. \W9(i : 
Mo. Man, Whitel . Il99fll : lOkoshi et aJ.1 . 120031 ^ sim 



ulations of galaxy formation at high-redshifts. 
All, these simulations indicate that DLAs are 
in majority faint (sub-L*) in small dark mat- 
ter halos Vc <C 100km s~^. Based on cross - 
section arguments, iFvnbo. Moller. fc Warred (^1999^ ■ 
Mealier et al.1 (|2002l 'l. and ISchavd (|200ll 'l arrived 
to the same conclusions. From the chemi- 



cal evolu tion point of view. JMatteucci et aL 
19971: I Jimenez. Bowen. &: Matteuccil . Il999l : 



Boissier. Peroux. &: PettinilJ2003r ) argued that DLAs 



are caused by gas rich low surface brightness 
dwarf galaxies, as seen locally in at least one case 
( Bowen. Tripp, fc Jenki^ l200lh . 

However, DLAs show asy mmetric profiles of 
their high ionization spe cies ( Prochaska &: Wolfd . 
19971 : iLedoux et al.l . 119981 ). This has been used 



to argue that DLAs ar e, in fact, d ue to thick 

Il98fil . Il 99,4 
others, e.g . 



massive rotating d i sks ( Wolfe et al 



Prochaska 



iyj2oool~ 



, -™- But 

I), iMcDo nald fc Miralda-Escud^ 



Mailer et aJ ■_ - 

(|1 9991 1 and |Haehnelt. Stein metz. Ranch I (|200nD . 

showed that a large range of morphologies can re- 
produce the observed kinematics: DLAs can arise 
from the combined effect of a massive central galaxy 
and a number of smaller satellites or filaments. In 
fact, cold gas accretion along filaments could be an 
important mechani sm, especially at high redshifts 



(KerisIeEaD,liQOi) 



Whether or not DLAs are indeed massive will 
lead to different clustering properties of the galax- 
ies around them. In hierarchical ga laxv formation 
models fe.g. lMo Rr. WhitFl . ll99fil . l2002 f). this clustering 
yields a measurement of the dark matter halo mass 
associated with DLAs relative to that of the galaxies 
used as tracers of the large scale structure. In partic- 
ular, if the galaxies are less (more) correlated with the 
DLAs than with themselves, this will imply that the 
halos of DLAs are less (more) massive than those of 
the g alaxies. Here, we use z ~ 3 LBGs (|Steidel et al. . 
Il999l ) as tracers of the large scale structure. 

In analyses similar to that presented here, 
Gawiser et al. I (|2no i1 found no clustering of galaxies 
aro und one single DLA at z = 4 towards BR 0951-04, 



and lAdelberger et al.l (|20o3) found a lack of galaxies 
near four DLAs (they found two within a cylinder of 
radius of 5.7/i~^Mpc and depth Wz < 0.025 whereas 
~ 6 were expected if the cross-correlation is the same 
as the galaxy auto-correlation). They argued that 



this is evidence that DLAs and LBGs "do not reside 
in the same parts of the universe" . It is important to 
note that both of these studies were not sensitive to 
scales larger than 5h~^ Mpc because because of the 
k^all field of view available. 

Other studies, however, have pointed to an over- 
density of galaxies near DLAs. IWolfe I 



(|1 9931 1 



com- 



bined several studies of Lya emitters around DLAs 
and found evidence for a correlation between emit- 
ters and DLAs a t a mean redshift z = 2 .6, signif- 
icant at > 99%. iFrancis fc HewettI ()l993l ) reported 



the discovery of super clustering of sub-DLAs at 
z ~ 2.4 and z ~ 2.9: a total of four Hi clouds are 
seen towards a QSO pair separated by 8', each being 
at the same velocity. Recent results from narrow- 
band imaging of the Francis & Hewett field shows 
that spectroscopically confirmed Lya emitters are 
clustered at the redshift of the strongest Hi cloud 
at z = 2.9 ( l ogA^H i = 20.9) towards Q 2138-4427 
( Fvnbo et all l200.1 'l. iBoche et al.1 (j200nl ) identified 
eight Lyman-alpha emitting galaxies near the DLA 
at z = 2.3 towards PHL 957 in addi tion to the previ- 
ously discovered Coup Fourre galaxv ILowenthal et al 



((199L 1. implying the presence of a group, filament, 
or proto-cluster associated with the DL A. Other ev- 
idence of cl ustering include the work of Ellison et al.l 
(l200lh and iH'Odorico. Petitieaii. Cristiarii {20o¥ . 
Ellison et aU ()200lh found that the DLA at Zabs = 
3.37 towards Q0201-I-1120 is part of a concentration 
of matter that includes at least four galaxies (in- 
cluding the DLA) over transverse scales greater than 
5/i-i Mpc. Id'Odorico. Petitiean. fc Cristianil (|2002l 'l 
showed that out of ten DLAs in QSO pairs, five are 
matching systems within 1000km s~^. They con- 
cluded that this result indicates a highly significant 
over-density of strong absorption systems over sepa- 
ration lengths from ~ 1 to 8h~^ Mpc. 

Despite numerous attempts to reproduce 
the pro perties of DLAs (e.g Katz et al 



199(il: iMcDonald MiraJda-Escude 



Mo. Mao, Whitel. Il 999l: iGardner et a1.l. 




Nagamine. Springel. fc Hernquistll2003bl:IOkoshi et al.1. 
2003h in numerical simulations or semi-analytical 
models, there are no predicti ons of the clu s tering 
of galaxies around DLAs. In iBouche et al.l (|2004l . 

in preparation), we discuss the cross-correlation of 
LBGs around DLAs in Smoothed Particle Hydrody- 
namical (SPH) numerical simulations. There, we find 
that the DLA-LBG cross-correlation is weaker than 
the LBG auto-correlation by approximately 30%, in- 
dicating that DLAs reside in less massive dark matter 
halo s than those of LBGs . 

In lBouche &: Lowenthall ( 20031 ). we presented a pre- 
liminary measurement of clustering of LBGs around 
a DLA at z ~ 3 towards APM 08279+5255, based 
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on part of this data set. We detected an over- 
density of S/Sg ~ 3 at the 95% level on scales 
2.5 < rg < 5 Mpc, implying that at least some 
DLAs reside in dense environments. In this paper, 
we expand on that study and present results on the 
DLA-LBG cross-correlation in our full survey of three 
fields. 

In sectionlJl we present the imaging data used here, 
our completeness limits and pre-selection of liigh- 
redshift galaxy candidates. We used photometric red- 
shifts to select our z = 3 LBG candidates as discussed 
in section 121 Our clustering analysis is shown in sec- 
tion|lJ The results on the DLA-LBG cross-correlation 
are shown in section \^ and discussed in section El 
Section [71 contains our conclusions. 

Throughout this paper, we adopt Q^j = 0.3, 
Qa = 0.7 and Ho = lOO/i km s~^ Mpc^^; thus, 
at z ~ 3, 1" corresponds to ^ 21.5/i^^ kpc and 1' 
to ^ 1.29h^^ Mpc, both co-moving. At that red- 
shift, H{z) ~ AAGHo, so 6z = 0.1 corresponds to 
67h-^ Mpc in CO- moving coordinates . 

2. THE DATA 

2.1. The DLA fields 

Given the allocated telescope night, we selected 
three fields for the presence of a DLA at z ~ 3 and 
with the additional constraint that the QSO must 
be at a higher redshift than that of the DLA, i.e. 
■Zabs << zqso- The redshift z ~ 3 is ideal for se- 
lecting LBGs efficiently using standard photometric 
filters. 

Our results on the QSO APM 0827 9 +525 5 field 
were presented in Bouche fc Lowenthall (|2003l ^ . The 
other QSOs in our survey are PC 1233+4752 
(zem = 4.447), with a DLA at ^abs = 
3.499 (|White. Kinnev. Beckel Il993l ): and SDSS 
J0124+0044 {zem = 3.840), with a DLA at Zabs = 
3.077 (C. Peroux, 2003, private communication). The 
individual QSO and DLA properties are listed in Ta- 
ble [H 

Although our DLAs do not meet the column den- 
sity threshold of logA'jjj > 20.3 often quoted, the 
latter is arbitrary and based on resolution thresh- 
old of previous-generation instruments. Furthermore, 
the metallicities and the Hi properties of 'sub-DLAs' 
are not different to the 'strict-DLA' population (e.g 



Peroux et al.l . 120031 ) . For the purpose of this study. 



the hydrogen column densities are logA'^jjj > 20.0, 
which ensure that the absorption is damped and that 
the gas is neutral. 

2.2. Observations 



The observations were carried out with the MO- 
SAIC camera at the Kitt Peak National Observatory 
4-m telescope on UT 2000 February 7 and 8 (run 
I), and on UT 2001 September 23-26 (run II). Run 
II was photometric, some cirrus were present during 
run I. The seeing was 0.9-1". 5 for both runs. 

The wide field imager MOSAIC contains eight 
2kx4k thinned SITe CCD. With 0".258 per pixel, it 
has a field of view of 35' on a side. The readout 
noise is ~ 6e~ pix~^, the dark current is negligible 
(~ 5e~ hr~^), and the average gain is 3e~ ADU^^. 
Each CCD has been thinned for detecting U-band 
photons. 

We imaged our three fields through four broad- 
band filters — U (Stromgren) and BV&I (Harris set) 
(see Fig. ^ — using a standard dither pattern (five 
pointings) to remove cosmic rays and detector de- 
fects. The total integration time for each field was 
typically 4hr(U), lhr(B & V) and 4hr(I); the obser- 
vations are summarized in Table |21 In addition, we 
observed several Landolt (1992) standard star fields 
through each filter. 



2.3. Data reduction 

The data were reduced with the package MS- 
CRED (v4.1) within IRAF3(v2.11.3), following the 
reductions guidelines of th e NOA O Deep Wide-Field 
Survey (Jannuzi & Greer], l2000l ) . This package was 
specifically designed to reduce MOSAIC data. The 
reduction process contains more steps than typical 
optical observations and is detailed here. 

We first performed the overscan level subtraction. 
We then corrected for the small correction (<0.5%) 
due to cross talk between adjacent chips. For each 
night's data, we removed an averaged zero frame, or 
bias frames, from the science images. The thinned 
MOSAIC chips required no dark correction since the 
dark current was only 5e~^ per hour. 

Flat-fielding is critical to achieve precise photome- 
try. However, we had to deal with two non-traditional 
complications: (i) the MOSAIC instrument at the 4- 
m suffers from a ghost image of the pupil in all bands 
due to reflections in the optics of the camera that 
needs to be removed; and (ii) dome flat-fields may 
match the night sky to only 1 or 2% (usually larger). 
Thus, sky-flatfielding in addition to dome-flatfielding 
was necessary and since the pupil image is an additive 
light effect, it has to be removed from the flat-fields 
first. 

In the dome flats, we removed the pupil image by 
fitting an axially symmetric pattern to the data them- 

: 0.268, r^A = 0.728, h = 0.71) changes these numbers by ~ 3%. 



Using the latest cosmological parameters from WMAP {Qm 
^IRAF is distributed by the National Optical Astronomical Observatories, which are operated by AURA, Inc. under contract 
to the NSF 
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selves with the task MSCPUPIL. We then flat-fielded 
all the science frames with pupil- free dome-flats. 

At this point, the pupil image was still present in 
the data. It is more difficult to remove the pupil 
from individual science images than from the dome- 
flats because the pupil pattern is much fainter, and 
the pupil image is mixed with all the faint and bright 
sources in the data. A simple threshold scheme to 
remove the objects is not feasible since the pupil is 
still present and the data poorly flat-flelded. Thus, 
to remove the pupil from the science images and to 
make the sky-flats, we had to extract the pupil image 
from the science images themselves through the fol- 
lowing iterative steps: (i) we created a sky-flat (ver. 
1) from the average of the science frames with a me- 
dian rejection; (ii) we extracted the pupil image from 
the sky-flat (ver. 1) using the task MSCPUPIL (pa- 
rameter 'type' set to 'data') and removed it from the 
science frames with RMPUPIL to produce a 'first- 
pass' pupil-free data; (iii) we created another sky-flat 
(ver. 2) using the pupil-free data and applied it to 
the data. We found that low-level light from bright 
stars was creating signal in the sky-flat even if strong 
min-max rejection was used. To solve this problem, 
we created object masks by using a 2g threshold on 
each of the eight CCDs and masked out large areas 
around the brightest objects. Steps (i) through (iii) 
were repeated using the masks and the flnal sky-flat 
was normalized and applied to all the frames. For the 
I-band onl y, we remo ved fringing using the procedure 
in (Jannu zi &: Gree"!! . I2OOO) before applying the final 
sky-flat. 

Cosmic ray removal was done using the task XZAP 
from the package DIMSUM and customized routines. 
Bad-pixel masks including the cosmic rays and bleed- 
ing regions were constructed. 

De-projecting the 8 CCDs on a single image is a 
two step process and requires very good astrometry. 
First, using the coordinates of several hundred USNO 
stars, we interactively derived astrometric solutions 
(RMS < 0.5") with MSCCMATCH for each dithered 
exposure. Then, we mapped the eight CCDs onto 
a single image by rebinning the pixels to a tangent- 
plane projection, thus producing pixels of constant 
angular size, with the task MSCIMAGE.'^ This pro- 
cess matches the World Coordinate Solution (WCS) 
solution of all bands to a common reference position. 

The individual flat-flelded, astrometrically cali- 
brated images with a uniform zero-point were aver- 
aged with an average sigma-clipping rejection to pro- 
duce the final stacked images. The scaling of each in- 
dividual dithered image was performed interactively 
on ~ 300 astrometric calibration stars common to all 

*The parameter 'fluxconserve' was set to 'no' because the science 
yield a constant sky per pixel. 



images. 

For the I-band of run II, we were not able to achieve 
a satisfactory sky-fiatfielding(residuals ~ 1%). In or- 
der to correct for this, we applied a median filtering 
to a block-averaged image of the stacked frame and 
applied the normalized results to the image. 

Even though all bands were deprojected to a 
tangent-plane solution using the same position and 
orientation on the sky, the relative pixel positions of 
objects in the different bands were not exactly iden- 
tical because of dithering and effects such as flexure 
of the telescope, and optical distortions due to fllters. 
Since it is important to have identical pixel positions 
for the photometry, we had to register and rotate 
slightly (~ 0.003 deg) each image with respect to the 
reference band (C/) using ~ 25 stellar objects with 
high S/N in all bands throughout the field of view. 

The rms in the relative astrometry of ~ 150 stellar 
objects with 20 < m/ < 22.5 and mu < 23.5mag is 
typically ~ 0.4 pix throughout the field. 

2.4. Calibration 

The standar d star frames, which contained ~ 150 
La,ndo]d (jl992l^ standard stars observed through each 
filter and airmasses 1.0 < X < 2.5, were reduced the 
same way, including the rebinning of the pixels to 
a tangent-plane projection. For run I, because only 
a handful of standard stars were available, we tied 
the photometric solution of run I against run II on 
the APM 08279+5255 field, which was re-observed in 
run II. 

Total flu xes of the stan dard stars were measured 
by fitting a lMoffatI (|l969l l profile (similar to the 7"- 
aperture flux used by Landolt 1992 at the 0.015 mag 
level). Table 01 presents the best fit (computed us- 
ing a linear Single Value Decomposition algorithm) 
to the photometric equation (e.g. Harris 1981): 



^obs 



c 



m"- = -2.5 log + ZP + aX + P (color), (1) 

texp 

where C is the number of counts, t^xp is the effective 
exposure time, ZP the photometric zeropoint, a the 
airmass coefficient and /3 the color term, which ac- 
counts for variations in the effective wavelength sam- 
pled by the filter for stars with different colors. We 
achieved a flux limit of typically (uiab) 25.5mag (3a) 
in a 2xFWIIM diameter aperture (see Table |1] for all 
the flux limits). 

2.5. Source detection 
All sources were detected in the I-band using 



the SExtr actor (v2.1.6) software (|Bertin fc Arnouta . 

exposures have been flat-fielded with flats equally distorted to 
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^9^. This algo rithm performs sour ce detection and 
photometry (see ISimard et all l2002l . for a summary 
of this package) . We optimized the configuration pa- 
rameters to ensure the faintest sources were detected 
and to optimize our completeness. We used the lo- 
cal background estimated in a 24-pixel wide annulus. 
The images were convolved with a 2 pixel FWHM 
Gaussian kernel before source detection. The de- 
tection threshold was set to 1.5 sigma with a min- 
imum area of 5 pixels. Bad-pixel masks are used as 
flag images. SExtractor is able to perform deblend- 
ing of close objects. The number of deblending sub- 
thresholds was set to 32 pixels, and through experi- 
mentation, the minimum contrast parameter was set 
to 0.0001. Our catalog contains approximately 40,000 
objects per field, 30,000 of which have / > 22.5mag. 

2.6. Photometry 

For each object, we measure the color in a 
2 X FWHM diameter aperture, where we took seeing 
variations on different bands into account: the color 
in two bands, e.g. {U - B) = mB(2xFWHMB) - 
m[/(2xFWHMt/), and similarly for other colors. Al- 
though this procedure is strictly valid only for star- 
like objects, it has been shown to be a g ood approxl 



mat io n for faint and unresolved galaxies (|Smail et al 
^9^. Indeed, from Hubble Space Telescope stud 



ies, the half-light radius of LBGs is, on average, 0.4" 
( Lowenthal et al. . ^93), much less than our seeing. 
For sources that were not detected in one band (i.e. 
flux< lo"), the magnitude in that band is set to the 
Icr flux limit in a 2 x FWHM diameter aperture and 
no color term is computed (Equation 1). 

In the remainder of this paper, we use aperture 
magnitudes. To convert those to total magnitudes, 
we estimate the total magnitude correction for star- 
like objects to be m/(tot) = m/ (2 x FWHM) - 0.35 in 
the I-band, calibrated by adding simulated stars with 
known total flux into our images and measuring the 
recovered flux in the chosen aperture. 

In addition, each object in our catalogs was cor- 
rected for Galactic extinction b y adopting E(B — V ) 



values taken from the maps of iSchlegel et al 
assuming an Ry = 3.1 extinction curve. 



2.7. Completeness 

In order to estimate our completeness, we added to 
our images fake stellar objects with MOFFAT profiles 
that matched the image point spread function (PSF). 
Fluxes were measured with the same aperture. We 
find that we are 50% complete up to mj ~ 24.35mag. 
The exact values for each field are shown in Tabled 
Using the transformation Jab = n^-i + 0.47, this cor- 
responds to Iab — 24.8mag (7?.ab~ 25mag, Steidel et 
al. 1993) and to 0.67L*, where ~ 24.5 for galaxies 



at z ~ 3 ( Steidel et al. . 19991 ). Our 90% completeness 
level is Jab — 24.4mag. Thus, we reached a depth 
sufficient to ensure that we sample well L* galaxies 
at z = 3. 

2.8. Selecting Lyman break galaxy candidates 

From our catalog of ~ 40, 000 objects, we rejected 
objects close to the field edges and objects with a 
FWHMJMAGE less than 0.85% the FWHM for stel- 
lar objects, presumably cosmic rays or bad pixels. We 
selected a subsample (1/3) using the following color 
cuts 

iU-B)AB > 0.3, (B - /)ab < 3.7, 
{U-B)ab > 0.3 + 0.55-[(S-/)ab-2], (2) 

and with 22.9 < /ab ^ 24.8mag. This removed most 
z < 1 objects and reduced the amount of computing 
for the photometric redshifts, which is explained in 
the next section. ~ 10, 000 objects per field met the 
criteria shown in Eq. |2 

2.9. Stellar contamination 

We did not perform any star-galaxy separation 
since beyond / ~ 23mag, number counts are dom- 
inated by the faint galaxies at those magnitudes. 
However, we determined empirically the stellar con- 
tamination at faint magnitudes the following way. 
We first select point sources, i.e. those with FWHM 
of the source < 1. 15 x FWHM of seeing, up to I < 
23.5mag and that fall in our color selection box 
(Eq. [21). We then select all the LEG candidates with 
22.5 < mj < 24.35mag that have a probability to 
be at z = 3 lb 0.25 greater than 50%, from our pho- 
tometric redshifts (see below). This subset contains 
some of the point sources between nij = 22.5mag and 
mj = 23.5mag, namely those with colors consistent 
with being at z ~ 3. Finally, we compare the num- 
ber counts of the LBG candidates with that of the 
point sources that have z ~ 3 colors and extrapolate 
the counts at mj > 23.5 to mj = 24.35mag, assum- 
ing the counts to be constant at m/ > 23.5mag (a 
conservative assumption given that the star counts al- 
most certainly decreases from mj = 22.5mag to mj = 
23.5mag). We found the stellar contamination to be 
< 7%, < 15% and < 13% for the APM 08279-^5255, 
the PC 1233-^4752, and the J0124-b0044 field, respec- 
tively. 

3. SAMPLE SELECTION & PROPERTIES 

Following our rough color-cuts, we used photomet- 
ric redshift techniques to select LBG candidates in 
narrow redshift slices. 
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3.1. Photometric Redshifts 

There are two approaches to photometric red- 
shift estimations: the e mpirical training set 



method (e.g. iKod . Il98,4 IConnollv et all . Il99 



and th e spectral energy distribution (S EP ) fit 
ting (iLanzetta. Yahil. Sz Fernandez-Sotol. 
Sawicki. Lin, fc YeA Il997l: I Budavari et al 



1996; 



200(1 : 



Fontana et al.l . l200ol : ICsabai et al.i . i200.¥ . The former 



is an empirical relationship between colors and red- 
shifts determined using a multi-parametric fit. The 
latter is based on a set of SED templates (empirical 
or theoretical). The two methods are comparable 
in their performance at z < 1; however, the train- 
ing set method is not always f easible, especially at 
high-redshifts (see discussion in Benited . 12000. ) . 



On the other hand, the SED fitting method works 
best when there is a strong feature in the SED, such 
as the 4000A break, or the 912A Lyman break. Thus, 
SED fitting methods were rap idly developed for the 
Hubble deep field (HDF) (e .g. iBudavari et"ZI . l2000l : 
iFernandez-Soto et al.l . l200lh with an accuracy of typ- 
ically Az ~ 0.06(1 + z). In our case, we used 
the co de Hyperz from Bolzonella. Miralles. &: Pelld 
( 200d ) . which includes intergalactic absorption due 
to the Lya forest and internal extinction Ay- We 
update d the intergalactic absor ption prescription fol- 
lowing iMassaro tti et al. ' ("2001) and we used the ex- 
tinction curve of iCalzetti et al.. (.2000. ) with Ay vary- 
ing from to 1.2. 

We used the templa te set made of the four 
empir ical SEDs of IColeman. Wu. WeedmanI 
(Il98d). extended in the UV (A < 1400A) by 
Bolzonella. Miralles. &: Pel ld (2000|) using the syn- 
thetic models of lBruzual^ Charloti (|l993l ) with pa- 
rameters (SFR and ag e) that matched the s pectra 
at z = 0. Note that iFernandez-Soto et alJ (|200ll ) 
exten ded the CWW te mplate s using the pow e r laws 



of .Kinnev et al 



^19931 ). and iBudavari et al.l (|200d ) 

used the extensions of iKinnev et al. I (Il99d ). The 

SED templates are convolved with the MOSAIC filter 
response curves (including the CCD response), and 
2phot is found from the maximum of the likelihood 
distribution C{z) derived from the distribution. 

3.2. Simple tests on the HDF 

We performed several tests of the technique on the 
sample of 150 spectroscopically confirmed galaxies at 
red shifts z < 6 in the Hubble deep field north (HDF - 
N) (|Cohen et al.l . l200d : iFW^andez-Soto et alJ . l2001^ . 
We experimented with 5 template sets that included 
the 4 CWW templates and various starburst tem- 
plates from Starburst99 ( Leitherer et al. Lll999l ). Of 
the 18 galaxies with 2.75 < Zphot < 4.5, two (11%) are 
outliers with Zspec — 1- From that sample, we found 
that (i) the four CWW templates gave the lowest 



scatter + Zspec) = 0. 053, both measured using 

the bi- weight estimator of iBeers. Flvn. &: Gebhardig 
(|l99d ). and (ii) A^/(l -t- •^spec) is not improved using 
near-IR photometry. Thus, near-IR observations are 
not required for photometric redshifts at z ~ 3. 

3.3. Using priors 

SED fitting methods give the most likely redshift 
given the observed set of colors. However, informa- 
tion such as size, or flux, can be included in pho- 
tometric redshift techniques using Bayesian proba- 
bilities (following .Benitez. . ,2000 ). We coupled the 
SED fitting scheme with the prior likelihood distri- 
bu tion for a. gala xy of magnitude mj parametrized 
bv lBenited (|200d ) in the following way: the product 
prior X likelihood is decomposed over the SED types 
T: 



P{z) 



^PT{z\mi) 
T 



Ct{z) 



(3) 



where pT{z\mi) is the prior probability given the 
galaxy magnitude m/, and Ct{z) is the probability of 
observing the galaxy colors if the galaxy is at redshift 
z and has a type T. The photometric redshift Zphot 
is taken from the maximum of the P{z) distribution, 
and the errors, u^, are computed from the FWHM of 
P{z) divided by 2.35. Redshifts with large Gz may 
be unrelia ljle. A good e stimator of reliability is the 
following (|Benitezl . E)Onh : 



Paz = P{\z - ZphotI < 0.2 X (1 + zphot)), (4) 

which estimates the 'goodness' of a photometric red- 
shift Zphot using Eq. |31 and also has the useful feature 
to pick hkely outliers (|Bemte d . l2000h . The factor 0.2 
is arbitrary, but since the rms of photometric red- 
shifts Gz is ~ 0.05(1 -|- z), this factor corresponds to 
approximately 4 x fj^. At z < 6, the overall rms of 
Az is 0.11, and Az/{l + Zs^cc) = 0.06, similar to 0.059 
found by Beni'tez (200Q.). 

3.4. Photometric redshift distributions 

Fig. 1^ shows the redshift distribution of the three 
fields. The dotted histogram shows the photometric 
redshift distribution using the CWW templates with 
no priors. The continuous histogram shows the pho- 
tometric redshift distribution using the priors. Fig. |21 
shows that using the priors has the effect of remov- 
ing galaxies with Zphot — 2 that are likely at lower 
redshifts. As expected, the distribution of galaxies at 
z ~ 3 is not affected much. This is due to the fact 
that this method is sensitive to the shape of the SED, 
which has a strong break between the U and B filters 
at that redshift. 
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We used the photometric redshift -Zphot of our galax- 
ies to determine their absolute magnitude M/^rest- 
The ii'-correction, which for galaxies at z ~ 3 cor- 
responds to the extrapolation of their intrinsic flux 
at Arest ~ 8OOOA from their observed flux at Arest ~ 
2OOOA, was computed using a weighted sum on each 
SED. Each template was weighted by the prior prob- 
ability pT{z\mf) • Ct{z) since the best-fitted SED was 
a combination of the spectral types T (see Eq. ISJ. At 
redshift 2; ~ 3, the ii'-correction is, however, small: 
it is typical lv ^ 0-2 for blue SEDs. su c h as the Irr 
template of IColeman. Wu. fc WeedmanI (llOSOT ). 

Fig. 121 shows the absolute magnitude Mj func- 
tion of z. Each dot represent one galaxy in our fields. 
The two continuous lines show our magnitude cuts 
and were computed using an Irr SED. The galaxies 
at redshift 2; ~ 3 are, as expected, between the two 
lines and near our completeness limit, which provides 
a consistency check of the photometric redshift tech- 
nique. 

Almost all points that are outside the range allowed 
by the continuous lines in Fig.|31are between the dot- 
ted lines, which represent the magnitude range for an 
E/SO template, and thus are best fitted by the E/SO 
type. This type has a strong break at 4000A that 
creates a large iT-correction of four magnitudes and 
makes these objects too luminous for their apparent 
magnitude. The fitted SEDs are likely to be wrong. 
At that redshift, Zphot ~ 2, IR photometry is needed 
to better constrain the SEDs. 

3.5. Selecting reliable galaxies in slices 

From the subsample described in section 12.81 we 
selected ~ 100 LBGs that are in a redshift slice cen- 
tered on the DLA. Specifically, they were chosen to 
have a high probability of being at the redshift of the 
DLA, Zs^hs, i-e. 



P(Zabs ± W,/2) = Pdla > 0.5, 



(5) 



where Wz is the redshift slice width. We also de- 
fined two additional redshift slices shifted by +0.15 or 
—0.15 from Zabsj P+ and P_, respectively. We choose 
a redshift width of W -,^ = 0.1 5 because, as discussed in 
iBouche fc Lowenthall (|2003 h. it produces the largest 
sample in the smallest redshift slice, given the rms 
of the photometric redshifts. At the end of § El we 
show that a different choice of Wz does not change 
the results. 

More importantly, this criterion (Eq. ^ corre- 
sponds to high quality photometric redshifts as illus- 
trated in Fig. H for the APM 08279+5255 field. The 
left panels of Fig.^lshow the probability distribution, 
and the right panels show the probability distribution 
as a function of Paz defined in Eq. ^ both for the 
three different redshift slices. Fig. 0fa), (b), and (c) 



show P_, Pqla and P+, respectively. The dots repre- 
sent galaxies detected in all four bands, UBV&: I. The 
filled squares indicate objects that are not detected 
in the U band. Smoothing the distributions using a 
Gaussian kernel (scaled to the peak) produced the 
continuous lines in Fig. ^ The dotted line shows the 
minimum threshold (> 0.5) used in selecting LEG 
candidates in each of the slices. From the right pan- 
els, galaxies that have a high probability of being 
in a redshift slice also have reliable photometric red- 
shifts, indicated by the fact that Paz > 0.9. For each 
redshift slice, the number of galaxies that met the 
threshold is shown in Table 13 

Fig. shows the x — y positions of our LEG candi- 
dates that met the criterion Eq. |S1 The square regions 
show the masks used to cover bright stars and defects 
such as streaks. 

4. CLUSTERING ANALYSIS 

In the curr ent hierarchical theory of galaxy forma- 
tion (e.g. see Longair 1 11998?. and reference therein), 
small quantum fluctuations that were stretched out to 
cosmological scales by inflation grew (mainly linearly) 
during the radiation-dominated era, till the present. 
The initial power spectrum {P{k) oc A:"), which char- 
acterizes these fluctuations in Fourier space, is nearly 
scale invariant (i.e. n ~ 1) on all scales. Initially, 
all scales grew at the same rate. Small scales en- 
tered the horizon before the universe became matter- 
dominated. During that time their growth was sup- 
pressed. The resulting power spectrum PE{k) has 
n ~ 1(— 3) on large (small) scales. These dark-matter 
fluctuations formed deep gravitational potentials in 
which galaxies and galaxy clusters formed. When 
the density contrast reached 5p/ p ~ 1, the fluctua- 
tions grew non-linear ly until A/9//9 ~ 200. 

Since only gravity is driving this build-up of matter, 
massive galaxies are more likely to be found in high- 
density regions, whereas low-mass galaxies are more 
uniformly distributed. This produces an enhance- 
ment of the clustering of massive galaxies. There- 
fore, the clustering properties of galaxies probe their 
dark-matter mass distribution. The auto-correlation 
^ (r) is a natural tool to study clustering in this con- 
text, since .^dm is the Fourier t ransform of the evolved 
power spectrum Psik) (e.g. iPeacock' 



I999I ). The 



galaxy correlation ^gg is related to the dark-matter 
auto-correlation ^dm via the bias b. At a given red- 
shift, 

^gg{r) = b'{M)^T,M{r). (6) 

This bias can be com puted in t he Press-Schechter for- 
malism extended bv IMo &: W hite (2002) and refer- 
ence therein. 

Similarly to the galaxy auto-correlation S^gg, one 
can define the cross-correlation ^dg between DLAs 
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and LBGs from the conditional probability of finding 
a galaxy in a volume d^ at a distance r = |r2 — ri|, 
given that there is a DLA at ri: 



P{LBG\DLA) = nu{l + idg{r)W2. 



(7) 



where Uu is the unconditional background galaxy den- 
sity. Thus, the number of neighbor galaxies in a cell 
of volume Ay is given by Np = N{1 + ^dg{''~))^ where 
= Ug/SV and ^ is the cross-correlation averaged 
over the volume Al/. This estimator of ^ requires an 
estimate of the unconditional background galaxy den- 
sity Ug. There are two ways to quantify Ug in case of 
the cross-correlation: one way i s to use gala xies spa- 
tially far from the DLA as in iBouche &: L owenthal 
(2003); the other way is to use the entire galaxy cat- 
alog. Here, we used the latter because of the simplic- 
ity of this method when randomizing the line of sight 
(see section 15. 2p . Naturally, large fields will yield a 
better estimate of Ug. 

We can extend the spatial cross-correlation ^ to 
angular correlation function w since the former is 
directly related to spatial correlation functions ^. 
Namely, if the selection function 0(2:) ~ 1/Wz within 
[—1^2/2,1^2/2] and zero otherwise, then 



w{r) 



fW./2 



(8) 



or w{rQ) ~ ^ (j^^ where /3 = 1 — 7 and A is 
a constant that can b e computed analy tically (e.g. 



a constant tnat can b e computea anaiy ti 
lAdelberger et al.l . 1200.11 : lEisensteinL |20o3). 

We used the following estimator of the cross- 
correlation Wdg{r), 



1 + Wdgir) 



Nr A'obs(r) 
A^gA^rand(r) 



(9) 



where Aobs(^) is the observed number of galaxies be- 
tween r — dr/2 and r+dr/2, A'rand('') is the number of 
randomly distributed galaxies, Ng is the total num- 
ber of galaxies, is the total number of random 
galaxies, and <> denotes the average over the num- 
ber of DLAs (A^dla)- The random catalog A^,- is much 
larger than the galaxy catalog so that the variance of 
A'rand(?^) is negligible. 

The errors to w i n Eq. IHl are given by (e.g. 
LauHvfc Sza]^.IT99^^ 



a,. 



Q"<A^obs> 

< A^obs > 



[1 + ^^^9], 



(10) 



where < Aobs > is the number of neighbors in the 
annuli averaged over A^dla- Using solely the Pois- 
son variance of the mean < A^bs >> i-e- c<Afobs> — 



\/ Anhs / A"r)T, A , yields the Poisson errors for Eg. El (e .g. 
Mo. Jing fc Boerneil Il992l : iLandv fc Szal^ . Il99.'^ : 



CT,, 



1 + Wdg 



V^DLA < A'obs > ' 



(11) 



In general, the cross-correlation and the galaxy auto- 
correlation will increase the variance to < Aobs > j 
and aw will be larger (see Eisenstein, 200^). 

Since < Aobs > is proportional to the total num- 
ber of galaxies, A^^, the expected rms of the cross- 
correlation function is 

1 



(T„, OC 



V^DLANg 



(12) 



as one might have expected. Thus, the noise in Wdg 
goes as the inverse of the square root of the number 
of DLAs, A^DLA; and as the inverse of the square root 
of the number of galaxies Ng. 

5. RESULTS 

With the clustering formalism layed out in sec- 
tion ^ we can present our results on the DLA-LBG 
cross-correlation (§ 15.1(1 and on the computation of 
the errors 15.2(1 . 

5.1. DLA-LBG cross-correlation from the combined 

fields 

Fig. El shows the DLA-LBG cross-correlation Wdg 
computed using Eq. |^ In computing Aobs('') 
and A^rand i^ Eq- ini we took into account the 
masked regions with bright stars shown in Fig. O 
The dotted line shows the auto-correlation Wgg of 



Adelberger et al.l ((2003) and the continuous line 
shows a fit to the amplitude of Wdg using Wgg as tem- 
plate. The fitting method is described below. 

It is necessary to take into account the different se- 
lections of the different fields in performing the sum 
in Eq.ini indicated by the brackets, so we must weight 
each field accordingly. We chose to weight each field 
according to its errors at each angular scale rj. Thus, 
for each field /, we compute -^^^ = ^+wi(r) in six an- 
null. Then, the combined angular cross-correlation, 
Wdg is computed from the weighted mean of wi{ri), 
where the weights were (Ti{ri) at each angular scale 

The error bars shown in Fig. o"(rj), are from the 
full covariance matrix to Wdg- The latter is found 
using the standard error propagation formula: 



GOV 



A^DLA 

E 

1=1 



df{n) dfjr, 
dwi{ri) dwi{rj 



-COYi{ri,rj) (13) 



where A^dla is the number of DLAs, /(r) is the 
weighted mean Wdg{r) and GOV; (rj,rj) is the covari- 
ance matrix of the individual field / given by Eq. 1141 
(below) . 
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5.2. Error computation 

Because each DLA is at a slightly different red- 
shift, each field has a different selection function, and 
it is necessary to take into account these differences 
by weighting each field accordingly. As mentioned, 
we chose to weight each field according to its errors 

The errors need to be computed carefully. Several 
options are available. The proper way to compute 
the errors would be to resample the DLAs (via boot- 
strap techniques), but this is impractical here given 
the number of DLA fields at our disposal. Another 
way would be to bootstrap the galaxies, which would 
reproduce only the Poissonian errors (Eg. 111(1 . 

We used yet another method, which is to perform 
Monte Carlo simulations in which we randomize the 
position of the DLA. This takes into account the clus- 
tering variance due to the galaxy auto-correlation, 
but misses the variance (and co-variance) due to the 
cross-correlation itself (the factor 1 + Wdg in Eq. HJ^ . 
However, this term will be small on scales larger than 
5h~^ Mpc because Wdg << 1. 

Thus, in each DLA field we computed the full co- 
variance matrix GOV; from = 200 randomizations 
of the DLA position: 



GOV 



1 



rir 



rir 



Wk{ri)-w{ri)][wk{rj)-w{rj)] 



CL4) 

where wt is the klh. measurement of w and w is 
the average of the measurements of the cross- 
correlation. The errors (Ji{ri) to wi{ri) for each field 
I follow. 

Our errors are consistent with the errors expected 
from our ana l ysis o f cosmological simulations: in 



Bonche et al.l (|2004l . in preparation), we conclude 
that with a data set of this size, we can be sensitive 
to the cross-correlation only on scales 5-10h~^ Mpc, 
which is where we see a positive cross-correlation. 

5.3. The integral constraint 

Because the unconditional galaxy density, in 
Eq. d is estimated from the total observed galaxy 
density, whereas it should always be lower than the 
observed galaxy density, all estimates of ^ (or w) are 
biased low. This bias Aw, often referred to as the 'in- 
tegral constraint ' , can be computed a nalyticallv (e.g. 
Landv Szd^ . 1199,4 ISa.sla.wL Eoooh . For the angu- 
lar cross-correlation function, it is : 



A-fw = ^y"dl7 Wdg{re), 



(15) 



where Wdg is a model of the cross-correlation and 17 is 
the total survey area. Thus, Aw will be smaller with 



larger fields, and we corrected Wdg in Fig. El for its 
bias C due to t he integral constra.int a, ssuming that 



Wdg IS Wgg 



from Adelberger et al. ( 20031 ) 



6. DISCUSSION 



We next discuss quantitatively how significant 
the cross-correlation shown in Fig. El is. We also 
show that our measurement of the DLA-LBG cross- 
correlation (i) is not reproduced by random lines of 
sight (at the 95% level) and (ii) is not seen in other 
redshift slices that do not contain the DLAs. We 
compare these results to past and future studies at 
the end of this section. 

6.1. Is this result consistent with no clustering? 

It is clear from Fig. El that we will not be able to 
constrain the slope of the cross-correlation, so we as- 
sumed that the cross-correlation w^^ has the same 
shape as the auto-correlation w^g. That is, we fitted 
the cross-correlation to the model 



(16) 



where a is the unknown amplitude. 

We fitted the amplitude of the cross-correlation a 
using the covariance GOVj computed in Eq. El us- 
ing the following x statistic: 



X 



E 



{re^ ) -Wdg {re, )] GOV, J [w^g {rg^ ) - Wdg {re, )] , 

(17) 

where w^g is the model of the angular cross- 
correlation and GOV~^ is the inverse of the covari- 
ance matrix computed using a Single Value Decompo- 
sition algorithm. The ^^{a) distribution is shown in 
the small panel in Fig. El In performing the inversion 
of the covariance matrix, we rejected the eigenvalue 
corresponding to the last radial bin or to scales sim- 
ilar to the size of the field where there is no signal 
(see discussion in Bernstein [1994]). 
We found that an amplitude of a > was favored: 



1.62 ± 1.32. 



(18) 



This measurement of the DLA-LBG cross- 
correlation is most significant on scales greater than 
5-10/i~^ Mpc. At this point, we cannot, however, 
conclude whether the DLA-LBG cross-correlation 
is stronger or weaker than the LBG-LBG auto- 
correlation. Taking — 4: h~^ Mpc and c:^ 0.6 
for the LBG auto-corre l ation ( Porciani &: Giavaliscd . 
'20021; lAdelberger et a,1.l . Hool), we find that the cor- 
relation length of the cross-correlation Wdg is Tq = 
5 ± 4.5/1-1 Mpc (or To = 7 ± 6.8/if/ Mpc). 
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6.2. Is this result drawn from random lines of sight? 

Given the large rms to the fitted amplitude a, could 
our result simply be a large fluctuation of the set of 
possible values for random lines of sight? To test 
this, we chose 100 lines of sight selected at random, 
excluding the central 5h~^ Mpc to ensure that the 
new lines of sight are not correlated with the ones 
centered on the DLAs. We then computed the cross- 
correlation for these 100 random lines of sight in the 
redshift slices centered on the DLAs. As before, we 
computed the weighted mean to Wdg and used Eq. 

Fig. [7| shows the logarithm of the x^(fl) fo'^ fixed 
amplitudes a for the 100 random lines of sight (filled 
circles). The continuous line shows the median of 
the distributions. The dotted and dashed lines are 
the 95% and 99% levels of the distributions. The 
median, 95% and 99% levels are found after a Gaus- 
sian kernel smoothing of the distributions using the 
optimum band width ( Wand 8z Jonesl 19951 ^ (The re- 
sults are not significantly changed using a fixed band 
width). The open square shows the location of the re- 
sult of Fig.ini Since it lies close to the 95% confidence 
level, this shows that the signal measured in Fig. El 
is not drawn from a random distribution of lines of 
sight, at the > 95% confidence. 

6.3. How about other redshift slices? 

The result of Fig. [S] should be compared with the 
cross-correlation when there is no DLA in the red- 
shift bin. From our photometric redshift analysis, 
we selected galaxies in two other redshift slices that 
did NOT contain the DLA. We chose the slices that 
were in the foreground and in the background from 
the DLA, and offset by -|- or —0.15 in redshift (see 
Fig. |1J). In each case, the fit does not favor any 
clustering: the best amplitude is a = —0.20 it 1.26 
and a = —0.24 it 2.04, respectively. A clustering sig- 
nal in this slice would have cast a strong doubt on 
our results that do contain the DLA in Fig. IHl In ad- 
dition, we performed the same check on another slice 
at redshift 3.6. The best amplitude for this slice is 
a = -0.13 ± 1.44. 

We repeated the analysis with Wz = 0.20 to 
test whether the observed clustering depends on the 
choice of the slice width. We found that a = 1.45 it 
1.35 in this case, so we conclude that the slice width 
does not strongly affect the clustering signal. 



6.4. Comparison with past and future work 
Wolfe I (|l993h also found that Ly-emitters are 



stronglv clustered around DLAs. In contrast 
Gawiser et al.l (|200ll ) did not find evide n ce of cluster- 



ing and the study of lAdelberger et al 
a lack of galaxies near their four 



( 2003t ) found 
DLAs, within 



5.7h^^ Mpc. Since these two surveys were not sensi- 
tive to clustering on scales larger than > 5h~^ Mpc, 
and ours is not sensitive to < 3 — 5/i~^ Mpc, our 
results are not inconsistent with theirs. The lack 
of galaxies on small scales could, however, be due 
to more local environmental effects, such as strong 
galactic winds from star forming galaxies. 
Altho ugh simula ti ons o f DLA properti e s exis t 

1^ iKatz et al.l (1199^): iGardner et a1.l ((2^3); 

iNa gamine. Springel. fc Hernquist' ('2003b')), no pre- 
diction of the D LA-LBG cross -c orrela tion has been 
published. In Bouche et al. (|2004l . in prepara- 
tion), we use the Tree-Smoothed Particle Hydro- 
dynamical (TreeSPH) cosmologic al simulations of 
iKatz. Weinberg, fc HernauistI (|l9961> ) to measure the 
'theoretical' DLA-LBG cross-correlation. These sim- 
ulations contain 128'^ dark matter particles and as 
many gas or star particles. Each galaxy (= > 64SPH 
bound particles) is able to form stars. With a sim- 
ilar number of DLA-LBG pairs and redshift depth 
(Wz = lll/i~^Mpc), we find Wdg > with the 
same signal to noise. Furthermore, with a much 
larger sample of 200 simulated DLAs, we find Wdg — 
0.75 ± 0.1, or ro ~ 3.5/^"^ Mpc. Thus, DLA halos 
are less n iassive than the halos of LBGs. which are 
lOifMe (|Porcia,Tii GiavaliscnI. Uml; lOuchi et al.. 
I2n03h . Given the present sample of three DLAs, our 
observed constraint on Wdg with its uncertainty is 
consistent with these simulation results. 

7. SUMMARY AND CONCLUSIONS 

Based on deep {fJ-i^ABisky) ~ 27.6 mag arcsec"^) 
wide-field images (0.31deg^ or ~ 65 x 65/iy/ Mpc 
co-moving at redshift z = 3) around three DLAs, 
we identify LBG candidates brighter than Iab = 
24.80mag using photometric redshift techniques that 
included the / magnitude as a prior estimate in ad- 
dition to the colors. 

From the redshift likelihood distributions, we se- 
lected LBG galaxies within a redshift slice of width 
Wz = 0.15(~ az) centered on the redshift of the 
DLAs 

^abs- Within that slice, we cross-correlated the 
LBGs with the position of the DLAs and found that 

• the amplitude of the DLA-LBG cross- 
correlation Wdg relative to the auto-correlation 
Wgg was Wdg/wgg = o = 1.62 it 1.32, corre- 
sponding to To = 5 lb 4.5/i^^ Mpc (co-moving), 

• the amplitude of the DLA-LBG cross- 
correlation is a > 0, which is significant at 
the > 95% confidence level based on Monte 
Carlo simulations, 

• the clustering signal was not present in three 
redshift slices that did not contain the DLAs. 
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Given the uncertainty of our results, we cannot 
put constraints on the halo masses of DLAs and 
discriminate betw een the large disk hypothesis (e.g 
Wolfe et all 198^^ . and small sub- L* hypot hesis (e^ 
Mailer et all 1200(11 : iHaehnclt. Stcinmetz. &: RauchJ 



20001 : iMc^her et al 



1^ 



2002) ■ Our observation of 
the clustering on large scales (> 4/i~^ Mpc) is 
not inconsistent with previous clustering s tudies 



( Gawiser et Zl . l200ll : lAdelberger et all 1200.1 ) 



since 



these were limited to small scales. In order to be able 
to directly compare these studies with our present 
results on scales < 4/i~^ Mpc, a larger sample of 
DLAs and multi-object spectroscopy of our LBG can- 



didates are needed. This will enable to test whether 
the cross-correlation is stronger or weaker than the 
auto-correlation. 
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Table 1 

Properties of the DLAs and the QSOs. 



APM 08279+525 PC 1233+4752 J0124+0044 



R.A. (J2000) . . 08''31'"41.6'' 12^35""31.r 0l''24™03.8^ 

Dec. (J2000) . . 52° 45' 17" 47° 36' 06" 00° 44' 33" 

/ 165° 45' 17" 130° 32' 12" 139° 58' 19" 

b 36° 14' 25" 69° 17' 27" -61° 02' 43" 

Au %mag) .... 0.20 0.08 0.13 

Eb-v OM 002 0.03 

QSO Properties 

QSO BAL 

ZQso 3.81 4.447 3.840 

niR 15. 2^ 20.63^ 17.9^ 

Radio >S'20cm = 1.3mJy ^ Seem < 77/iJy ^ 520cm = O.llmJy 

DLA properties 

2abs 2:974 3A99^ 3.077' 

Hi Wr (A) . . > 4.8^ 4.22^ 

logTVjji (cm"^) < 20.36 20.9^ 20.1^ 

[Fe/H] -2.31*^ 



^Galactic extinction from iSchlegel et al.l (|1998l ). averaged over the field. 



(l2002l^ (^) Sclineider. Schmidt. & Gunn 


' ("1991') (A) 'CariUi et al.1 (12001) 


(5) IWhite. Kinnev. & Becker 


(ll993^: (6) .Petitiean et al.. (.20001^: (7) Peroiix. 



C., 2003, private communication. 



Table 2 

Summary of the observations. 





U band 


B band 


V band 


I band 








Field 




Total Exposure Time 




UT Date of Obs. 


APM 08279+5255 


3.75hr 


35min 


50min 


2.08hr 


Feb. 


7, i 


i, 2000 


PC 1233+4752 


3.50hr 


40min 


50min 


1.92hr 


Feb. 


7, i 


i, 2000 


J0124+0044 


3.72hr 


47min 


52min 


2.08hr 


Sept. 23-26, 2001 
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Table 3 
Photometric solution. 





Run I 




Run II 




Filter 


ZP 


ZP 


a 




U 


23.26 (0.08) 


23.52 (0.02) 


-0.421 (0.02) 


0.018 (0.007) 


B 


25.085 (0.03) 


25.26 (0.02) 


-0.180 (0.01) 


0.095 (0.008) 


V 


25.07 (0.02) 


25.25 (0.02) 


-0.072 (0.01) 


-0.025 (0.009) 


I 


24.58 (0.05) 


24.76 (0.02) 


0.020 (0.01) 


0.003 (0.011) 



^For run I, we assumed the airmass coefficient a and the color term 
P to be the same as for run II. 



Table 4 
Depth of the observations 







Exp. /Frames 


Airmass 


FWHM 








Fields 


Filter 


(sec./#) 


(min-max) 


(arcsec) 


(mag/m^s) 


(mag/m^s) 


(mag/m^fi) 


APM 08279-F5255 


U 


13500/15 


1.07-1.17 


1.1 


27.78/28.49 


26.03/26.74 


25.85/26.56 




B 


2100/7 


1.14-1.27 


1.1 


28.44/28.36 


26.69/26.61 


26.49/26.42 




V 


3000/10 


1.21-1.37 


1.2 


28.21/28.23 


26.46/26.48 


26.17/26.19 




I 


7590/20 


1.11-1.62 


1.1 


27.15/27.61 


25.40/25.86 


25.21/25.66 


PC 1233-F4752 


U 


12600/14 


1.04-1.16 


1.05 


27.82/28.52 


26.07/26.78 


25.94/26.65 




B 


2400/8 


1.05-1.07 


1.0 


28.59/28.51 


26.84/26.76 


26.75/26.68 




V 


7590/10 


1.07-1.11 


0.9 


28.18/28.20 


26.43/26.45 


26.45/26.47 




I 


6900/15 


1.07-1.45 


1.1 


27.19/27.64 


25.44/25.90 


25.30/25.76 


J0124-F0044 


U 


13400/16 


1.19-1.51 


1.5 


27.88/28.59 


26.13/26.84 


25.60/26.31 




B 


2800/7 


1.24-1.42 


1.5 


28.69/28.61 


26.94/26.87 


26.44/26.37 




V 


3100/7 


1.46-1.90 


1.4 


28.33/28.35 


26.58/26.60 


26.14/26.16 




I 


7500/20 


1.21-1.80 


1.1 


27.67/28.13 


25.93/26.38 


25.75/26.21 



^The airmass is where ( is the zenith angle of the telescope. 

''Limiting surface brightness in magnitudes per square arcsecond. 
'^Measured inside a 2x FWHM diameter aperture. 



Table 5 

Number of galaxies in the different redshift slices. 



Field 


P_ 


Pdla 


P+ 


APM 08279+5255 


89 


84 


17 


PC 1233-F4752 


20 


70 


83 


J0124-F0044 


184 


65 
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Fig. 1. The solid lines show the transmission curves for our four filters U, B, V, and I. The dashed line shows the CCD 
response function. The dotted lines show the filter transmission convolved with the CCD response function. 
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PC 1233+4752 




phot 



APM 08279+5255 




J0124+0044 




>hot 



'phot 



Fig. 2. — Redshift distribution for each of our fields. The dotted histogram shows the photometric redshift distribution using 
no priors and the template set A. The continuous histogram shows the photometric redshift distribution using the priors. Using 
the priors has the effect of eliminating the large number of galaxies that have been assigned Zphot — 2 wrongly, but does not affect 
the distribution at ~ 3 significantly. The vertical dashed line shows the redshift of the DLA zula- This plot shows the effect of 
the priors and that our selection peaks at a redshift close to that of the DLA. 
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PC 1233+4752 





Fig. 3. — Absolute magnitude vs. photometric redshift for each field. Each dot represents one galaxy in our fields. The 
rest- frame absolute magnitude Mi was computed using the distance modulus and the it'-correction (sec text for details). The 
continuous lines show our magnitude selection 22 < 7ab < 24.8mag and are computed using an Irr SED. At z ~ 3, galaxies are 
between the two continuous lines and near our completeness limit as expected, which gives us more confidence in the photometric 
redshifts. The dotted lines show our magnitude selection for an E/SO SED. Clearly, points that are outside the range allowed by 
the continuous lines are best-fitted by the E/SO type, which has a strong break at 4000A and thus a large iiT-correction. This plot 
shows that the photometric redshift technique is self-consistent. 
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Fig. 4.— Redshift slices centered (a) in front of, (b) on, and (c) behind the DLA for the APM 08279+5255 field. The value 
^abs=2.974 is indicated by the vertical dashed line. Each dot represent a galaxy that was detected in the four UBV& I bands. 
The filled squares indicate objects that are not detected in the U band. The left column shows the probability distribution as a 
function of photometric redshift. The continuous line shows the smoothed distribution (arbitrarily scaled to the peak). The right 
column shows the probability to be in that particular slice as a function of the 'goodness' of the photometric redshift Paz defined 
in Eq.0] The dotted line shows the minimum threshold (50%) used in selecting LBG candidates in each of the slices. 
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Fig. 5. — For our three fields, the x y position of our LEG candidates relative to the QSO location. North is left, East is down. 
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Fig. 6. — The cross-correlation Wdg between DLAs and Lyman break galaxies in a redshift slice of width {Wz ~ 0.15) that 
contains the DLAs. The filled squares show the cross-correlation for the combined fields. The dotted line is the LBG auto- 
correlation Wgg (from AdclberEcr ct al., 2003, using Eq.|H]to account for the volume of our redshift slice). The continuous line is 
a fit to the amplitude of the cross- correlation using Wdg = a x Wgg, i.e. we assume that both Wgg and Wdg have the same slope. 
The small panel shows the distribution as a function of the amplitude a and the la range. 
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Fig. 7. — The points show the logarithm of X^i^-) a function of the amplitude a for 100 random lines of sight (excluding the 
central 5h~^ Mpc) in the redshift slice centered on the DLAs. The continuous line shows the median of the distributions. The 
dotted line and the dashed line are the 95% and 99% confidence levels, respectively. The open square shows the location of the 
fitted amplitude shown in Fig. El which shows that the signal measured in Fig.|^is not drawn from a random distribution of lines 
of sight at the > 95% confidence level. 



